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DEPENDENCE  OF  REACTIVITY  OF  THE  WATER  BOILER  ON  THE  MASS  OF  25*  IN  THE  SPHERE 


By  J.  Hinton 


ABSTRACT 

A  semitheoretical  determination  of  the  percentage  change  in  reactivity  of  the  water  boiler  per 
gram  of  25  removed  from  the  sphere  is  presented.  The  method  consists  in  calculating  the  external 
neutron  source  required  to  keep  level  the  neutron  intensity  in  the  boiler  when  x  grams  of  25  are  re¬ 
moved  from  the  sphere.  Use  is  made  of  the  experimentally  determined  neutron -distribution  curves.  It 
is  found  that  18.2  g  of  25  subtracted  from  the  sphere  produces  1  per  cent  change  in  K. 

*  *  ♦  *  * 

The  evaluation  of  many  experimental  results  obtained  with  the  water  boiler  requires  a  knowledge 
of  the  change  in  reproductive  factor  per  gram  of  25  removed  from  the  sphere. 

Several  experimental  methods  of  determining  the  dependence  of  reactivity  on  grams  of  25  have 
been  suggested  by  certain  members  of  the  theoretical  department  in  order  to  avoid  the  trouble  of  calcu¬ 
lating  the  dependence  theoretically.  One  of  these,  the  bpron-bubble  method,  was  carried  out  by  Kerst 
and  de  Hoffman.  It  seemed  worth  while  as  a  check  on  this  method  to  carry  out  also  a  semitheoretical 
determination.  In  the  calculation  here  presented,  use  is  made  of  the  experimentally  observed  neutron- 
distribution  curves.  This  makes  the  procedure  much  simpler  and  perhaps  somewhat  more  reliable 
than  a  complete  theoretical  determination. 

If  K  is  defined  as  the  number  of  fissions  produced  per  generation  by  each  fission  in  the  boiler, 
then  if  K  equals  1,  each  fission  produces  on  the  average  one  more  fission,  and  the  boiler  runs  at  con¬ 
stant  power. 

When  the  boiler  is  set  to  run  with  K  equals  1,  the  quantity  to  be  determined  is  the  change  in  K  if 
X  grams  of  25  is  removed  from  the  sphere. 

If  for  a  moment  the  spacial  distribution  of  the  activity  is  disregarded,  the  boiler  may  be  run  at  a 
constant  rate  even  when  K  is  less  than  1  by  adding  a  source.  If  this  source  is  equivalent  to  ^  fissions 
per  imit  time,  then  the  total  number  of  fissions  F  will  be  given  by 

F=:0/(1-K) 


or 

5K  =  l“K=:(p/F  (1) 

This  formula  can  be  used  for  calculating  the  change  in  K  due  to  removal  of  x  grams  of  25  from  the 
sphere.  For  this  formula,  first  it  must  be  determined  what  source  must  be  added  in  order  to  keep  the 
intensity  level  whenx  grams  is  removed  so  that  the  system  becomes  undercritical.  Since  both 
fissions  F  and  source  ^  are  actually  functions  of  r,  proper  averages  must  be  taken  in  calculating  6K. 
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In  order  to  determine  the  added  source  necessary  when  x  grams  is  removed,  assume  that  the 
boiler  is  critical  with  m  grams  of  25  in  the  sphere;  then  the  number  of  fissions  per  cubic  centimeters 
would  be 

F^(r)  =  n(r)  v(m/V)  (2) 

where  Of  is  the  fission  cross  section  per  gram,  and  V  is  the  volume  of  the  sphere. 

If  X  grams  is  removed,  leaving  the  volume  unchanged,  the  neutron  flux  as  well  as  the  density  of 
25  will  change.  Assuming  we  start  with  the  same  number  of  fissions  in  the  first  generation,  we  shall 
then  get  in  the  second  generation. 


FgW  =  [n(r)  -  e(r)]  [(m  -  x)/V]  va^  =  Fj^(r)  -<^(r) 


(3) 


where  €(r)  represents  a  small  variation  in  neutron  density. 

In  order  to  keep  the  intensity  unchanged,  namely,  to  have  F2  =  Fj,  an  external  source  has  to  be 
added,  distributed  in  such  a  way  as  to  produce  the  fission  distribution  (/>(r).  If  F(r)  and0(r)  were  every¬ 
where  proportional,  then  5K  would  simply  be  equal  to</>(r)/F(r).  Actually  it  results  that  0  and  F  have  a 
very  different  dependence  on  r.  Nevertheless,  a  formula  similar  to  Eq.  1  can  be  applied  as  follows. 

A  source  of  fission  neutrons  is  put  in  a  boiler  that  is  almost  critical.  The  boiler  acquires  a  steady 
intensity  that  depends  on  the  position  at  which  the  source  is  put  and  which  has  a  maximum  when  the 
source  is  at  the  center.  Now  defined  is  a  function  0(r)  equal  to  the  ratio  of  the  intensity  when  the 
source  is  placed  at  a  distance  r  from  the  center  to  the  intensity  when  the  same  source  is  placed  at  the 
center. 

Thus  by  multiplying  any  source  placed  at  a  position  r  by  the  function  0(r),  the  source  is  determined 
which  it  would  be  necessary  to  place  in  the  center  in  order  to  make  the  boiler  run  at  the  same  intensity. 

For  a  source  of  Ra-Be  neutrons,  the  function  0(r)  was  determined  experimentally  by  Kerst.  His 
results  are  plotted  in  Fig.  1.  It  is  to  be  expected  that  the  function  0(r)  is  somewhat  dependent  on 
energy  of  the  neutrons  emitted  by  the  source  and  is  flatter  for  neutrons  of  high  energy.  For  fission 
neutrons  having  a  lower  average  energy  than  the  Ra-Be  neutrons,  a  somewhat  steeper  curve  is  ex¬ 
pected.  On  the  other  hand,  it  is  expected  that  the  fission -neutron  curve  would  be  slightly  flatter  than 
a  thermal-neutron  curve.  The  desired  curve  would  therefore  lie  somewhere  between  the  Ra-Be  curve 
and  a  thermal-source  curve  of  the  function  0(r).  It  can  be  shown  that  the  thermal-distribution  curve 
is  very  nearly  equivalent  to  a  thermal -source  curve.  It  is  exactly  equivalent  with  an  infinite  water 
tamper. 

The  calculation  has  been  carried  out  both  by  taking  6^(r)  =  <5(r)Ra-Be  ^  where  n(r)  is 

the  experimental  thermal-distribution  curve.  The  result  is  practically  independent  of  the  choice. 

As  stated  above,  a  source  placed  anywhere  can  be  reduced  to  an  equivalent  source  placed  at  the 
center  by  multiplying  it  by  0(r).  Consequently,  the  source  distribution  <j){r)  will  be  equivalent  to  a 

2 

concentrated  source  of  strength  Air  I  0(r)r  0{r)dr  at  the  center.  And  in  a  similar  way,  the  source 

•^0 
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distribution  F(r)  is  equivalent  to  4Tr  F(r)r  0(r)dr.  Equation  1  can  now  be  used  m  the  form 

/•R 

/  <^(r)r20(r)dr 

6K  = - - - 

r  F(r)r20(r)dr 


The  fission-source  distribution  F(r)  is  assumed  to  be  proportional  to  the  thermal-neutron  distri¬ 
bution  n(r)  according  to  Eq.  2.  It  is  therefore  obtained  from  the  experimental  thermal  curve  (Fig.  1). 

The  remaining  part  of  the  calculation  consists  in  determining  the  function  <^(r).  Since  according 
to  Eq.  3,  this  depends  on  €(r),  the  change  in  neutron  density  as  a  function  of  radius,  €(r)  must  be 
determined  first. 
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An  equation  in  6(r)  may  be  obtained  by  taking  the  difference  of  the  diffusion  equations  with  m  and 
(m-x)  grams  of  25  in  the  sphere. 

It  is  assumed  that  the  concentration  of  hydrogen  is  practically  imchanged.  Therefore  the  scatter¬ 
ing  mean  free  path  and  the  slowing-down  process,  both  of  which  depend  mainly  on  hydrogen,  are 
assumed  not  to  change  when  x  grams  are  subtracted. 

Since  only  the  absorption  mean  free  path  and  the  neutron  density  change,  the  following  equations 
may  be  written,  where  index  s  refers  to  “sphere,"'  and  index  t  refers  to  “tamper” 

With  m  grams  n  =  n(r)  1/A  1/ A 

s  s 

With  m-x  grams  n  =  n(r)  +  €(r)  1/A  -  1/A^  -  n 

s  s 

Sphere  Boundary  Tamper 

With  m  grams  (X  v/3)An-  nv/A  +  q  =  0  *  (X  v/3)A  -  nv/A  +  q^  =  0  (5) 

s  s  s  L  n  V  V 

With  m-x  grams  (XgV/3)  A(n”  €)-(n-~  e)(l/A^~7?)  +  q^  =  0  ||  (X^v/3)A{n  ~  e) (n  -  €)v/A^  +  q^  =  0  (6) 

(  V32-in.  stainless -steel  sphere  wall) 

With  the  boundary  condition  that  the  number  of  neutrons  diffusing  through  the  stainless-steel  wall  from 
the  sphere  is  equal  to  the  number  diffusing  back  into  the  sphere  from  the  tamper  plus  the  number  ab¬ 
sorbed  in  the  wall 


% 
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where  is  the  absorption  cross  section  per  square  centimeter  of  the  wall  (0.015  cm^/cm^  for  a  y32- 
in.  steel  wall).  Subtracting  Eq.  6  from  Eq.  5  the  equation  in  (r)  is 

Ae(r) -[e(r)/L^j  =  n(r)  T)  ||  A€(r) -[eW/L^]  =  0 

2 

where  L  =  AX/3.  The  boundary  condition  for  6(r)  is  the  same  as  that  for  n(r). 

Substituting  U(r)  =  e(r)r  X/3?7,  the  equations  to  be  solved  are 


U"  (r)  -  =  rn(r) 

l2 


U"  (r)  - 


0 


U'(R+)  =  ^  U'(R-)  +  U(R) 


Xt  R 


(7) 

(8) 

(9) 


U(0)  =  0 


(10) 


-4. 


Taking  Lj  =  30  cm,  =  2.69,  =  .5,  =  1800,  1/Ag  =  0.075  +  10  x,  and  Xg/S  =  0.153,  the  solution  is 


-ar 

rn(r)e”^^  dr  -  — 


rn(r)e+^’^ 


dr  -  1.4914{e^*^  - 


(11) 


where  r  <  R,  and  a  =  1/Lg.  The  integrals  in  this  expression  are  evaluated  numerically,  taking  the  experi¬ 
mental  Mn- distribution  curves  as  n(r).  U(r)  being  known,  0(r)  is  calculated  by  Eq.  3.  Then  6K  is  ob¬ 
tained  by  numerical  integration  of  Eq.  4. 
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By  this  method  it  was  found  that  18.2  g  subtracted  from  the  sphere  produces  1  per  cent  change  in 
K.  The  result  of  the  boron-bubble  experiment,  as  given  in  LADC-816,  is  in  exact  agreement  with  this 
calculation. 


RADIAL  POSITION  OF  200- MC  SOURCE,  IN. 

Experimental  thermal  curve. 
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